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I. INTRODUCTION 

The description of the motion of extended bodies in 
Einstein's theory of gravitation is a complicated and 
many-faceted problem. Nearly all applications of Gen- 
eral Relativity crucially depend on our ability to describe 
how matter moves under the influence of the gravitational 
field. 

When it comes to the description of extended bod- 
ies, one usually has to resort to the use of approxima- 
tion schemes due to the complexity of the theory. In 
this work we utilize a multipolar approximation method, 
originally devised by Tulczyjew [1|, to characterize the 
motion of extended test bodies. We explicitly work out 
the equations of motion at the monopolar, dipolar, as 
well as quadrupolar order with the help of this method. 
In doing so, we put particular emphasis on the defini- 
tion of multipole moments, as well as on the definition 
of combined quantities, e.g. the mass or the spin of ob- 
jects, at different orders. The process of canonicalization 
is carried out in detail, and the resulting equations of 
motion are worked out in a fashion which allows for an 
easy comparison of quantities at different orders. 

Without going into historical detail we would like to 
point out, that multipolar methods were among the first 
methods to be studied in the context of the problem 
of motion in General Relativity. In particular Tulczy- 
jew's method can be traced back to the seminal work of 
Mathisson Q. We only note in passing, that several of 
the pioneering concepts of Mathisson's approach to the 
problem also resurface in other (later) works. For a more 
detailed account on the history of different multipolar 
approximation schemes we refer the reader to [3j. 



Although we solely focus on the description of the mo- 
tion of test bodies, the results obtained here are also rel- 
evant for other approximation methods. Most important 
arc the post-Minkowskian (3411 and the post-Newtonian 
approximations - see, e.g., the reviews [7|-|9( and refer- 
ences therein. The former method is useful to describe 
the scattering of an unbound and the latter method is 
commonly used to describe the inspiral of a bound sys- 
tem of two compact objects. For both approximation 
schemes the test mass results in the present paper can 
be extended to self-gravitating objects, if one relies on 
certain regularization techniques. Furthermore, our re- 
sults can be used as input for higher order black hole 
perturbation schemes [l(HlH or other radiation reaction 
calculations [131 1 14] . All of the above mentioned meth- 
ods are used to model gravitational waves emitted from 
different astrophysical sources and are therefore of direct 
relevance for the field of gravitational wave astronomy 
0. 

The structure of the paper is as follows. In section 
HI1 we briefly recapitulate the basic ingredients of Tul- 
czyjew's multipolar scheme. This is followed by the de- 
compositions of a general set of moments in section Mil 
which are crucial for the derivations in the succeeding 
sections. In sections IIVI - I VII we explicitly derive the 
monopolar, dipolar, as well as quadrupolar equations of 
motion. In IVI Bl a detailed comparison to the multipolar 
approximation scheme by Dixon is performed. We draw 
our final conclusion in section ["VTT1 Appendices [AI IB1 and 
[C] contain a brief overview of different quantities and our 
conventions as well as some useful transformation rules. 



II. BASIC DEFINITIONS AND THEOREMS 
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Conceptually the multipolar method of Tulczyjew [l| 
is based on the assumption, that the motion of an ex- 
tended test body along a representative worldline can be 
characterized by a set of multipolar moments, which are 
built from the energy-momentum tensor T ab of the body. 
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As in the case of other multipolar approximation 
schemes in the context of General Relativity, the start- 
ing point is the covariant conservation of the energy- 
momentum tensor, i.e. 



V h T ab = 0. 



(1) 



The general strategy consists of working out the con- 
straints of this equation on a general set of multipolar 
moments, which are covariantly defined via an expansion 
of the energy-momentum density of the following form: 

/+oo 
{t ab s {4) + v c [t cab s {4} ] 
-oo 

+V d V c [t dca % ) ]+...}ds. (2) 

Here the t ab are general multipole moments, (5(4) = 
<5(4) (x a — Y a ) is the delta function, Y a (s) characterizes 
a representative worldline of the body, and s denotes the 
proper time along this line. In other words, the con- 
tinuous energy-momentum density T ab of the extended 
body under consideration is replaced by an infinite set 
of moments, defined along a worldline - which is com- 
pletely arbitrary at the current stage. It is natural to 
expect, that the general energy-momentum conservation 
law (fTJ), imposes constraints on the moments defined via 
(|2|) , which in turn characterize certain properties of the 
body and its motion through spacetime. 

Of course, there is no simplification of the general prob- 
lem of motion at the present stage. Insertion of the full 
representation © into JT]) just yields a very complicated 
differential equation for the moments t abcd ---, which is ac- 
tually of infinite order. The approximation comes from 
the fact, that one truncates the series in @, and only 
considers a finite number of moments. The general ex- 
pectation is, that certain features of the body and its 
motion are adequately characterized by a small number 
of moments. This is of course analogous to multipolar 
expansion techniques as used in the context of classical 
mechanics or electrodynamics. In the following, the or- 
der of approximation is characterized by the notion of 
single-pole, dipole, or quadrupole, depending on which 
of the moments in the expansion @ are retained. 

On a technical level the method of Tulczyjew [l| is 
based on a generalization of the Du Bois-Reymond the- 
orem - called theorem B in the following, cf. section 3 
in [l|. Before stating theorem B we need to define what 
Tulczyjew calls the canonical form. An arbitrary singu- 
lar tensor density A Cl '" c " is said to be in canonical form 
if it can be written as 



k=Q 



[a Cl - Cfcbl - 6 "<5 (4 ) (x a -Y a )] , 



where the coefficients a Cl 



c 1 ...c k b l ...b n 



C!...C k b 



U c , a 



satisfy 

(c 1 ...c k )b 1 ...b n 
0. 



a 

ci...e fc &i...fc 



(3) 
(4) 



Here we used the abbreviation u a := dY a /ds for the 
tangent vector along the worldline. As was shown in [l6| 
it is always possible to achieve the canonical form. 

Theorem B states that, if for a tensor density A x " n 
and an arbitrary tensor field T^...^ we have 



A bl - b "T, 



bi...b n 



= 



(•5) 



in an arbitrary four-dimensional region D, then all the 
coefficients a Cl '" Ch ' 1 " n of the canonical form of the den- 
sity A bl --- bn vanish. 

For further reading on the method of Tulczyjew, we 
refer readers to lecture V in Hal. 



III. DECOMPOSITIONS WITH RESPECT TO 
THE VELOCITY 

Let us consider a set of general moments t abc - - with 
the following symmetry properties: 



j.c-L...c n ab £Ci...c n (ab) 

±C\...Cn<xb ^{c\...c n )ab 



(6) 



Note that the symmetry in the second line is motivated 
by the integral representation of moments in the multi- 
polar approximation scheme of Papapetrou [l7j , see also 
[18| . In the context of Tulczyjew's formalism there is 
a priori no reason to impose the symmetry in the first 
n indices. One could carry out the calculation without 
imposing it. However, this would only lead to an unnec- 
essary complication of the derivation, for the contribu- 
tions from the antisymmetric parts of t Cl - Cnab could be 
absorbed in the course of the canonicalization procedure. 
We provide an explicit example of the absorption process 
in section IVH 

With the help of the projector p% := 5% — u a Ub we may 
decompose the moments from ^ with respect to the four 
velocity, remember u a u a = 1, in the following way: 



t ab 


= °o ab + 2>uV 


°o« 


■= t cd p a c u d , 


°ab 


— t cd n a n b 



t 


:= t cd u c u d . 



(7) 

For the dipole moment, or three index quantity, we get 

t abc = o abc + 2o a ( fc U C )+O a uV + M a t 6c , 

> := t de fp a d u eUf , 
l ab := t de fp a d p b u f , 



■= t P d P e Pf, 

\bc +dbc 

t .— t Ud- 



(8) 
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For the quadrupole, or four index quantity, we get 

t abcd = labcd + 2 ab(c U d) +l ab U C U d 



-u a u b 2 t cd + 2u {a t b)cd , 



j^ab +efqh a b 

o := t JJ p e p f u g u h , 



o 



abc iej : gh a b c 



1 PePfPg u h, 



'abed ,ej ' gh a b c 



it.. / un a u c „a 
t PePfPgPh, 



t cd ._ t efcd UeUf ^ 

\ acd ±eacd 

l . - — L LLq . 



(9) 



The decompositions in ([I])-© are going to play a central 
role in the upcoming derivations. 



IV. SINGLE-POLE 



If we consider a single-pole object we start from 

o = v a f ab = /v„ [*<%)]■ (10) 



In order to bring this equation into the canonical form, 
we decompose the integrand with respect to the first in- 
dex, i.e. the one which is contracted with the derivative. 
Here and in the following we abbreviate the orthogonal 
projection of an index with respect to the velocity by 
a hat over the corresponding index, e.g., t abc 
The decomposed version of (|10[) then reads 



p b t adc_ 







V a [(t & » + u a u c t cb ) 5. 



(4) J 



This integral can be split up by means of 



V a [u a T c 



as follows: 

= 



ds 



\u r t 



cbl 



J (4)J 



5(4) 



ds 



5(4), 



(11) 



(12) 



(13) 



This form of the integral allows for the application of 
theorem B, i.e. the equations of motion are now given by 



— [u c t cb ] = 0, t hb = 0. 



(14) 



If we insert the orthogonal decomposition of t ab as given 
in (|7|) - note that in the single-pole case we could have 
started right away with this decomposition, without mak- 
ing the intermediate step in (|11[) - the equations of mo- 
tion in (fl4l take the form 



S_ 

ds 



o a + u a t 



= 0, 



o ab + o a u b = 0. 



(15) 



From the second equation, due to the orthogonality, we 
can infer that 



o ab = and o a = 0, 



(16) 



which leads to - after reinsertion into the first equation 
in OB 



o S 

t= const and — u" = 0. 

ds 



(17) 



In other words, we have shown that the equations of mo- 
tion of a single-pole particle take the form of the geodesic 
equation. Equation (|17[) suggests to identify the quan- 

o o 
tity t with the mass m of the test body. This result 

is of course not new, with the method outlined above, 

it was already derived by Tulczyjew in [l[. Finally, we 

note that the corresponding singular energy-momentum 

tensor is then given by 



T 



ah 



n b c 

m u u 0(4). 



POLE-DIPOLE 



(18) 



At the pole-dipole order we start from 

= V„f ab = h a [t ab S {i) ] + f V Q V C [t cab S (A) ] . (19) 



In order to be able to apply theorem B, we need to bring 
(fT9"f into canonical form. First, we focus on the second 
term with the two covariant derivatives. Following the 
procedure outlined in the preceding section on the single- 
pole particle, we start with a decomposition of the indices 
which are contracted with the derivatives, i.e. 



V a V c [t cab S {4) ] = / V a V c [(i £ab + t Mb u a u d ) 5, 



(4) J 



+ / V Q V C [{t d ' ab u c u d + u c u a u d u e t deb ) 6, 



= I ^(t deb Udu e )S (i) 



5u a 



u d u e t deb 6, 



ds— " (4) 

4b 



Va 

V a V c (t m ,S 



(4)J 



(4) 



(4)J 



+ / V a V c (t cdb u a u d 5 {i) ) 



(20) 



The last term in this equation can be rewritten with the 
help of 



Vy-7 rjiabc 1 / ) crjiabd 
[aVfeji = -Kabd 1 , 



(21) 



and subsequent application of ([12)) as follows: 
V a V c [t £db u a u d ] S (i) 



= v 



(4) 



+ / R ace b u a u d t cde S {4) .(22) 
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In order to arrive at the canonical form, one still needs 
to symmetrize the indices appearing in the second last 
term of (f20|) . This again produces a curvature term at 
the lower order. Combining the rewritten form of (|20[) 
with the results at the single-pole order from the previous 
section, equation (|19[) takes the form 



(2t £de u a u d + t £he ) 



'(4) 



Su a 



ds 



(j-dab , j-adb \ , " "* j.. 
(t Ud + t Ud) + -j^-UdU e t 



deb 



+t &b ] %)}+ / V a V 



t (ca)br 



(4) 



= 0. 



(23) 



This is almost the canonical form, we still need to ensure 
the orthogonality of the first two terms in the second last 
line in ([23]). Once again we apply (|12l) and obtain 



ds 



(t Mh u d + t &dh u d ) 5 



(4) 



5(4) 



5u c 
ds V ds 



u d t (cd)b ) <5(4). 



(24) 



Reinsertion into (|23p finally yields the canonical form of 
(HT 



cdb 



u c u d 



ds 



t Cb U r 



1 

2 Ra 



(2t U Ud + t ) 



2-—u d t 

ds 



Hi) 



(cd)b 



+ <a 1 <%)} + / v - v 



(4) 



= 0. 



(25) 



With the help of theorem B - starting at the highest 
order - we can infer from (1251) that 



(26) 



= t {£&)b =o< ca ) b + o< CQ V. 



Here we made use of the decomposition © in the last 
step. Transvection of (|2"B1) with the projector pf and the 
velocity Ub yields two conditions, i.e. 



o (ca) =0, and o< e 



0. 



(27) 



From the last equation, together with the identity 

Icab = 1 {ca)b + l Q {bc)a_ ^(ab)c^ ^8) 

we can infer that 



At the second highest order theorem B yields 
S (.<„a\u \ 5u a 



2p °Ts( t{cd)bud ) ~ is- UcUdtdcl ,+ 



t ab = 0. 



(30) 



Insertion of the decompositions from (|7J) and © leads to 



ds 



o a u b = 0. 



(31) 



Multiplication by Ub and reinsertion of the result yields 
two equations. These allow us to express parts of the 
orthogonal decomposition of the single-pole moment in 
terms of the decompositions of the dipole moment as fol- 
lows: 



5 



-UdP c ~T~ \° + u + t 



\c n ,d i ,cd 



o" u = -PdPcj; I + u 

Taking the antisymmetric part of (|33[) yields 
We introduce the spin in the following way: 



(32) 
(33) 

(34) 
(35) 



Note that the prefactor is conventional, in particular the 
minus sign comes into play because we started with a 
positive sign in front of the dipole term in (p~9|) . Now 
Q34j) turns into the well-known equation of motion for 
the spin @,[l7j], i-e. 



a b SS cd 
PcPd — 



S S ab 
ds 



u u c - 



6 S cb 
ds 



l 

Sj qac 

u b u c ^- = 0.(36) 



Furthermore, if we make use of the first equation in (|27[) 
and the definition of the spin (|35p . we can express parts 
of the orthogonal decomposition of the dipole moment in 
terms of the spin and the velocity as follows: 



- S ab u b , 

1 1 u 1 r n 

-~S ab -u c S c[a u b l 



(37) 
(38) 



From the lowest order in (|25l) we get, again via theorem 
B and by insertion of the decompositions from ([7]) and 
©, the following equation: 1 



, ,u d —t db + o b + tu b -u c o 
as \ ds 



cb 



\c„,b 



iff 6 

2 -fi-ace 



2u 



a (o ce +o c u e ) + o cae +o' 



0. (39) 



; cab 



o. 



(29) 



1 Here we introduced the shortcut 



ds 



5 



Taking into account the symmetries of the quantities in 
this equation and our findings in ([27]), @, <E3), 
([37]) . and (|38l) . we can rewrite (|39j) as follows: 



— p" 

ds 



0. 



(40) 



This is the equation of motion for some kind of general- 
ized momentum, which we define by 



p" 



• 1 c 

t -u c u d S 
8 



to u 



ud-r- S 

ds 



u c u d ^- t cd ) u b 
ds 



h,i 



u d -r- S 

ds 



(41) 



The second line serves as a definition of the mass to which 
now - in contrast to the result at the single-pole order - 
contains also contributions from the spin as well as from 
the transversal component of the decomposition in (|8|). 

Our equations of motion for the pole-dipole test body 
in this section are the most general ones. No a priori re- 
strictions were imposed on the decompositions in (|7J and 
(0. Furthermore, it should be stressed that no assump- 
tions were made regarding a possible spin supplementary 
condition. 

Equations (1551) and (|4"0j) are nowadays usually called 
the Mathisson-Papapetrou equations. In particular the 
characteristic spin-curvature coupling at the dipole-order 
is already present in Mathisson's pioneering work 
Note that Mathisson's equivalent to equation (14TJ1) has 
a slightly different form. This is due to the fact that he 

sets the o a component in the orthogonal decomposition 
of the dipole moment to zero - in other words he makes 
use of a supplementary condition - at an early stage in 
his calculation. On the other hand, Papapetrou does not 
impose any supplementary condition in his derivation in 
[l 71 ] . The equations of motion given by him are formally 
equivalent to (j3"o) and fiOj) , but his moments are defined 
in a different way, cf. [3] for more details. 

Finally, let us derive the energy-momentum tensor at 
the pole-dipole order. In terms of the spin, we have 



h cd +h c u d +t cd = -\s 

hence (j3"2")) and (|53"j) become 



cd 



+ u e S e{c u a >+ t 



i cd 



UdP c 



8 1 cd 5 1 
-r S cd --rt 
ds ds 



cd 



i 8 1 j 
b a ,cd 

-p d Pc Ts t 



pi a u b) S c 



+ ^ dP a c S ca , (42) 
'u e . (43) 



If we use this result - as well as all the constraints on 
the components of the decompositions in (J7J) and JS]) ob- 
tained in this section - in (|19[) . with the help of (fT"2|) the 
singular energy-momentum tensor for pole-dipole parti- 
cles becomes 



rpab 



L [a p b) S (4) - 



S c{a u b} 6, 



(4) 



(44) 



Note that (|4~4")) is not in canonical form. 



Supplementary conditions and conserved 
quantities 



The system of equations in (|36|) and (j40|) is under- 
determined. This is evident from the appearance of the 
projectors in equation (|36|) . Thus, supplementary con- 
ditions, or constitutive relations, involving the spin are 
needed to close the system. Before we discuss the im- 
pact of different conditions, we rewrite the equations of 
motion as follows: 



1 -a 

m u 



S i 
— m 
ds 



1 ij 8 { Si 

e ri dc d ci a t n 

-u S tidce -Pb-r\ u c-rS 

2 ds \ ds 



-u r — l u b S 
ds 



be 



— s ab = 2pI Q it b l. 
ds 



'"' 1 ,(45) 
(46) 
(47) 



The first two equations are obtained from the orthogonal 
decomposition of (|4T)j) . 

There are basically two covariant supplementary con- 
ditions at the pole-dipole order which have been studied 
in the literature, i.e. 



S ab u b = (*), 



S ab k = 0. (**) (48) 



To our knowledge, the first condition can be traced back 
to an early work of Frenkel and the idea for the 
second condition appeared first in a work by S y n g e (20j 
in a special-relativistic context, see also 13. [2l| - [24J ] . For 
both conditions there exist constant quantities, namely 



ds 



l 

rn 



ds 



5 i 

— TO 

ds 
S 



6 

ds 



PaP a 



S i i i. 

__ rt oab 

— , O ab o 

ds 



= o, 



= 0, 



(49) 
(50) 
(51) 



Note that without the imposition of any supplementary 
condition the derivative of the alternative mass parame- 
ter m fulfills 



1 1 S 1 6 P a 1 S 1 h 

m m—m=—^-P b — S ab - 
ds ds ds 



(52) 



Furthermore, in case the background spacetime allows 
for a Killing vector field ip a , the quantity 



S i 8 /i li,, 

Ts e:= Ts r ^ + 2 s v ^ 



(53) 



is conserved. For other (non-linear) conserved quanti- 
ties at the pole-dipole order see [H, [2(| and references 
therein. 
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VI. POLE-DIPOLE-QUADRUPOLE 

At the pole-dipole-quadrupole order the variational 
equation takes the form 

= V a f ab = J V a [t ab S {i) ] + J V Q V C [t cab 6 {i) ] 

+ J V a V d V c [t dcab S (i) ] . (54) 

In order to bring this equation to canonical form, we fo- 
cus on the third term in (|54|) and proceed along the same 
lines as in the single-pole as well as in the pole-dipolc 
case. With the help of the projector, the quadrupole 

2 Note that one has to be careful when it comes to the usage of 
the " " notation in combination with the hat " " notation for 



moment can be decomposed as follows: 

t dcab = t d£&b + u a t dceb Ue+u c u d t gfab UfUg 

+u a u c u d t 9feb u e u f u g 

+u a (u d t f£eb + u c t dfeb )u e u f 

+ (u d t e££lb + u c t de£lb )u e . (55) 

Due to their length, we provide the canonical form for 
the separate terms in (|55p. As in the previous cases, 
the canonical form is achieved by repeated application 
of (fT2|) . and the generalized version of (f2Tj) for multiple 
derivatives. 2 

projected indices. 



J V a V d V c [^ ab 5 (4) ] = J V a \7 d \7 c [t^ b S {i) }+ J \7 c ^R ade b t £dde S {4) + ^R ade £ t £ddb 5 {4) 



l 5 r 

3 ds 



Rade 



c_, +eddb 



5(4) + ■^Rdae b ;ct Cdae 8( 4 ) > . 



(56) 



V a ^d^c[u a t dCeb U e S {4) ] = I V Crf V c) \ p C eP d f — 



5 r 



t feab l 



5(4) \ + / V J 2R acf b u a t d£ ^u e S {4) + R acf d u a tf' ceb u e 5 



(4) 



-2p a e — [u c t ecaD u a \ 5 (4) } 4 



_5_ r 

ds 



Racf d U a U d tf £eb " 



5(4) " Racf b ;dU a t £def U e S {i) 



+2 



ds 



u d u c t d£ab u a 



5(4) " Rdce b U d Uft^ ae U a S {4) \, 



(57) 



/ 



V a V d V c [u c u d t 9fdb u f u g S {4} ] = J V (a V d) {u d t< 3 ^ b u f u g 8 {A) } + J V«{p^ [t 9f " b u f u g ] 5 



5(4) + -RadeVt 9 ^ 



/ "S°(4) f > 



(58) 



J W a V d V c [u a u c u d t^ e 



W a \ 3ii a — [t 3feb u e u f u g ] 8 (i) +2p a h u h t^ eb u e u } u g 5 (i) 
[t^ eb u e u f ug\ 5(4) - 2— [u a u a t 3feb u e ufu g ] 5(4) 



ds 3 



+R adc b u d u a t^ ec u e u f u g S (4) 



(59) 



V a V d V c [u a {u d t^ eb + u c t d f eb )u e u f 5 {i) ] = J V (c V a) {2u a tf £eb u e u f S (4) } 

+ J V c |2p»A_ [t^ UeU/ ] 5(4) - R gd /u a u d t^ eb u e u f 5 (i) 



2— (E oc /u a )^ ees We u /0(4) 



UeUfSn) + 2 — 

as 



; rfs 2 ^ 



''(4) 



+3R acg b u a — [t f£e9 u e u f ] 5(4) + R cd g b . a u a u d tf £e9 u e u f S {i} L (60) 



7 



V a V d V c [(u d t ecab + u c t deab )u e S (4) ] = 



V( d V c) |2p^- 

d„Jj.aceb 



aefb 



»(4) 



V d <j R fce a un aceb u a 5 {i) + R fce b u } t acae u a 5 {i) - 

u c udt a{&<i)h 



ds 

6„,d„\ +acfe 



'(4) 



a ( c 



(5 (4) - 2R dce b u d ii f t ac ^u a 5 



(4) 



3(4) 



(61) 



Together with our result in (|2"5j) . equations (|51))) - (j6"Tj) 
yield the canonical form of (I5H) - which we do not dis- 
play here explicitly due to its length. As in the previous 
sections, the next step in the derivation of the equations 
of motion consists in the insertion of the orthogonal de- 
composition from ([T])-©- To save us some work at the 
quadrupole order, we are not going to use the general 
decompositions in (JZ])-© directly in the following, but 
rather transform them to a more compact form first. The 
new decompositions shall be given by 



n ab + 2n( a uV+nu a u b , 



(62) 



t cab = ^cab + 2 ^c(a u b) + ^c u a u b 



i_dcab 



[dcab + 2n d < a u b) +n dc u a u b . 



(63) 
(64) 



As becomes apparent from (|6"3")l and (|64p . the 
terms in © and © have been absorbed in the 
composition. The explicit transformation laws 
the old and the new decomposition read 



1 2 
t and t 

new de- 

between 




n 



n' 



t +U e Uf 



r 5t ef _ S 2 t ef 2 _5_ f St 9ef 
ds ds 2 ds\ 9 ds 



2u 9 t cd{e R gc d f) 



o a + P>f 



o ab + p a eP ) 



1 c 


1 c 


P>S U 9 ( 


n 


= o c + 


li ca 


= o ca -| 


- PlP^g ( 


n cab 


= l cab 


+ P»\ 


n ab 


= l ab , 




2 cab 


= l Cab , 




^ dcab 


^dcab 



r sl ef _ s 2 t ef + 2 _s_ f St gef 

ds ds 2 ds\ 9 ds 



2 f 

St efa 



ds 



2u a t cd{e R gcd f) 



(65) 
(66) 
(67) 
(68) 
(69) 

(70) 

(71) 
(72) 
(73) 



r 



Hence, the energy-momentum tensor is given by 



T ab = / t ab § 



(4) 



V c [t cab 6 (i} ] + / V d V c [t dcab 6 {4) ] 



(74) 



From the form of the decompositions of the moments in 
(l62"T) - (f6"4)) and in (|55)l it becomes apparent, that only a 

small fraction of n ooc " terms contributes to the canonical 
form of VbT ab , when we make use of the re-defined t abc --- 
moments. In particular, only the integrals in (|56p and 
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(I57p yield non- vanishing contributions. 



The canonical form of the derivative of (l74l) becomes 



V a V d V c {t^ a %)}+ f v d v 



(c d) 
PePf- 



ds 



t-feab 



Jcd)b 



'(4) 



^Ts ^ ahUa ~ 2ilct " abua ) + Rac f b { tdm + 2w Q t d " £e/ u e ) + R ac / Q t /£afc + u a t^ e! 1 



^db 



'(4) 



A ^i? ac / w a Mrf t^ eh Ue + ^R ade c u c t iddb + 2u d u c t d£ab u a + t cb u c - u c u d t cdb 

-\Rdae h ^~ ae ~ Raef" ; d U a t^ U e ~ R dce "u d U ft^ U a + ]-R ac b (2t Me U a U d + t^) 



5(4) y = o. 



(75) 



Application of theorem B and insertion of the decompo- 
sition from (|64|) , yields at the highest order: 



2(dca)6 + l(dca) u b = q 



(76) 



Orthogonal decomposition with respect to the open index 
b leads to two symmetry relations, i.e. 



^{dea)b = Q and &(<*«*J=Q. 

The first relation in (l77l) allows us to infer that 



^dcab -^abdc 



(77) 



(78) 



Equation ([5H)) can be used to rewrite h abc as follows: 



'ds 



= ^Ts ( 2 

ds 



fcgf _ lfe u g_ ^gf u e + ^eg^ 



bd u c + lcad u b + ^cbaA ^ 



; dc-b i bd 



cb-cl 



(82) 



Note that in the second step we made use of the second 
relation in (|77p. which yields an algebraic symmetry for 



n abc , i.e. 



abc , 2 bca+ z cafc = ^ 



hence n abdc has symmetries similar to the ones of the 
curvature tensor 3 . 

At the second highest^ order, again by application of An analogous relation holds for « abed 
theorem B to (|75p. and insertion of the decompositions 
from (|63|) and (|64|) . we obtain 



(83) 



(" /6b+ " /e " b ) + " (Cd)b+ " (Cd)ub = °- ( 79 ) 

The orthogonal split of (|75t yields two "constraint equa- 
tions" - here and in the following we are going to use 
this name for equations, which allow us to express cer- 
tain parts of moments in terms of parts of higher order 



moments, e.g. n a 



, (cd)b 



be 



i abc (lobe ^ab 



namely: 



„d b 



-PePjUb 



l(cd) 



n / e s+ n fe u 3 



(80) 
(81) 



Equation (|81[) can be simplified as follows: 
2„rw . S 2 



fic: 



(84) 



Let us now turn to the next order, i.e. the second line in 
75]) . Applying the same procedure as before we obtain: 



\db , "d,.b 



+R.<t'pt 



1 /2 



as 



n e0 + h e u° - 2 n ec0 + rj e V it 



n fcab + n fca u b I + I n 1 L "+ n 1 "u" I « 



3 /<*-L W« 



2 rfca / + 2 dca u / , 2 <fc/ + ^c u /\ a 



= 0. 

(85) 



The orthogonal split of (|85p yields two equations 



2 2 

3 Using 177> . and remembering n dcab =n (dc)(ad) ^ one can check 

that n adci, + n cadf, + n Mca + n cMa = -2 n dcab holds. On the 

other hand, one also has n dabc + n bdac + n cabd + n bcad 
abdc _ >p ne j e £^ nan( j sides f these equations are identical, thus 
([78} holds. 



2 

2 H 



n db 



~d b 



S ri 



-p-p- ^ + A V - 2(n ec ^+ n ec u f )u c 
-Racf 9 \p b g n dcaf + S b n dca u f + 2p b l dcf u a 



+26 b n dc u f u a + - nf cab p d + n f cb p d u a 



(86) 
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n = - Pe u b - 



fn ecb +n ee u b> j 



n eb +n e u b -2[n ecb +n ec u b )u 



1 



/ca 



(87) 



The antisymmetric part of (|86[) can be viewed as the 
direct generalization of (f34|) . and suggests that the spin 
at the quadrupole order should be defined as 



a b 



s 



yielding 



1 [ab] + ^[a u b] _ 2 (^c[afa] + ^[a^j^ 



= -\pip)^r S ef + Racf 9 \pf n d ^f + 4" 



+2pl 6 n*^" + 2S [b n d ^ c u 1 u a + I n^pf) 
y y 3 



- h^ b pfu a 
-ln d b — hf - [d o b] R f £ 

2 f ds Pg^acf 



+4n ea f'u c + 2n ea u f u c 



4 2 
— n 
3 



eacf 



(89) 



(90) 



In contrast to our findings at the dipole order (f31)|). this 
propagation equation for the spin contains - as expected 
- also contributions from the quadrupole moment, which 
couple to the curvature of spacetime. 
From (EH) and (1511) we obtain 



-Rac 



a (h ce + h c u e ^ + i (h cae + h ca u e ^ 



u f u a (nf ce + rV c u e ) 



= 0. 



(93) 



This equation is the analogue to (j3"9"f found at the pole- 
dipole order. With the help of (|82 |) . (|87l) . (|9il. and (|92l) 

, 2 

it can be rewritten in terms of S and the n a — from 
the decomposition of the quadrupole moment in (|64p . In 
order to keep the equations at a manageable size, we 
introduce the auxiliary quantity 



A ab := h ab +h a u b -2(n cab +n ca u b ^ju c (94) 

2 i 
■ S al 



2 as 



with the properties 

u a A ab 
u a A ab 

u b A ab 

j^[ab] 

u a u b ^-A ab 
as 



= 0, 

= -u a A ab , 

= n a - 2 n CQ ?i c =u 6 5 
1 2 >, 

rr ao 

- 2 S , 

= S ah u a u b . 



(95) 



(96) 
(97) 



5 a 



(99) 
(100) 



This allows us to rewrite (l86l) and (l87l) as follows 



n a = 2n ab u b -S ab u b , 



(91) 



: ab 



$ 2. 



1 2 



as 2 

-2(n ' 



(n c6a +n c ( a u fc ))u c , (92) 

which can be viewed as the generalizations of (|37j) and 
(|38|) to the quadrupole order. 

What remains to be analyzed is the lowest order in 
([75)) . Application of theorem B at this order yields 



ds 

+ 1 -u d ( 

-u c (n cb + h 



- 2u d u c (n dc 



n ecab + n eca u b 



•+ n dc u b 



+R ace :d 



n'-u ) + n u + n 

2 (2 j 2 



b\ 

n u > 



- (n dac 
3 V 



+ n dac u e ) - u 



a ^ 2 cdc 



+ n cd u' 



db 



_ A (db) + u (b A d)e 

A< h u d) 

-A e fu e u f u b u d + pi d p b jR acf 3 

- n fcae + 2 n eac u f + 2 n ea u f u c 
3 



(101) 



ds 



- l fcae + 2 n eac u f + 2 n ea u f u c 
3 



+2u<> d p b Juf 



5A e f 5A^ 



ds 



-p d e u b —A eb - R acf 9 
ds 



ds 



n dcaf u n 



+2 l dcf u a u g + i n fca p d + lf c p d u a 
3 



Hence equation (1931) turns into 



(102) 



(103) 



ds 



nu b +n b -A cb i 



Racd He 



'l 2 

— n 

.3 



dcab 



i 2 C 
- n 

3 



; deb 



u a + n ac u"u" 
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- fi> b 

n b.x 



*>ace 



(I 


^ dcae 


2 2, 

f - n dca 
3 









i^e u a + ^«fc u e u oj + ^^^(2 n ce u a + 2 k c u e u a + A cae + n ca u e ) 



(104) 



c/.s 



n u b - p b e u d A ed - A cb u c + R ace d u d ( ~ n ecab + i n eca u b + n ec V+ n ec w b u' 

- ibce,.a.. i .2 eca _6 , 5 ec „&„.a 



-i?ace d ( ™ 6cae "d + 2 n bce u a W(i + - n eca p b d + n ec p b d u 



l -R 1 

2 1 L ace 



2w a u [e S c]f u f - u a S ce + u e u a S cf u f - 4u a ii d (n dec + n d< - c u e h - 2u a ^- n< 

ds 



1 2 

^u e S ca 
2 



+2u e it Q fu d S dc + 2u d n dc J + 2— n aec + 2 (n dec u a + n adc u e + n aed u c ^ u d - n ca ii e - n ec ii a + n ae u c 



+2u e u d n dca - 2 (nf ce + ^ c u e ) u a u f 
S(A bd u d ) 



_ p 6 



n w b + 2A^u d - ■ h //, </,///' 



,<L4 e 



+i?ace rf ( | n abce « d + 1 n ec Vu d + 4 n bae u c u d + n ec u a u b u d - ~ n ec V b - n ec u a p b d 



2 2 dcae + 2 2 d ca u e + ^dce u a + fadc^a 

3 3 



1 2 



(105) 



ds 



ace 


2 






[(•"• 


+ 




2 


a ds 


5 



1 2 



(5(n eca ) 5(n ec u a ) 



ds 



ds 



D 6 

^ace \d 



2 2 dcae 2 2 dca 2 dc 2 ^ 

3 3 



n + 5 ac u a w c - -R ace d n eca u d u 



p b e qclc td b 



/ 4 2 9 4 2 2 > 

> i d d \ ZL abce , ,1 oae c , ^ aec b { c\ TLae c b 

2 2 rfcae + 4 2 dce ^ a + 4 2 aed u c + n° 6 uV 

3 3 3 



(106) 



(107) 



This concludes the derivation of the equations of motion 
at the quadrupole order. Equation (|107p replaces ([3T)| as 
the new center-of-mass equation of motion. 

Analogously to the pole-dipole order, we introduce 

2 

combined quantities for the mass m and the generalized 

2 

momentum P a at the quadrupole order, as follows: 

m := n + S ac u a u c -^R ace d n eca u d , (108) 



P b := ^ u b + Ua A|ba + i?ac /g l^ Ud 

+4 n bae u c u d + | n aec p b + 2 n ae ti c p^ (.109) 

With these definitions the propagation equations for the 
spin and the center-of-mass (|107[) take the form 



c 2 

n a n b — Q cd — n^n^Rj A 
ds ~ 3 



- n edcf + 8 n edf u c + 4 n ed u f u c 
3 



J 2, 
P £ 

ds 



2 2 

— n 
3 



dcae 



4 2 7 4 2 

, dce„,a I r, aea,,c 



i ae n ,c n ,d 



3 3 



(110) 



(111) 



Our equations of motion for the pole-dipole-quadrupole 
test body in this section are the most general ones. It 
should be stressed that no assumptions were made re- 
garding a possible spin supplementary condition. 

Before we introduce a new combined quantity for the 
quadrupole moment, we make contact with the multipole 



formalism of Dixon in section IVI B[ This will allow us to 
bring the equations of motion (|110p and (jllip into a very 
compact form. 

From the system of equations (|110l) and (|lllj) it be- 
comes clear, that the evolution of quadrupole compo- 
nents is not constrained - in the sense that there is no 
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dedicated propagation equation for the n ah ~. Equation 
(PJ seems to generate only equations of the constraint 
type at higher orders. Furthermore, one is already forced 
to introduce a supplementary condition at the pole-dipolc 
approximation in order to obtain a closed system of equa- 
tions. Hence, it is rather natural to expect, that addi- 
tional supplementary conditions, now also involving the 

higher order moments n ah ", are needed at the quadrupo- 
lar order. The choice of such conditions depends on the 
type of body under consideration. We are not going to 
touch upon the question of possible choices for such a 
condition in this work. 



A. "Non" -conserved quantities 



In this section we calculate the derivatives of the 
22 2 
masses m and m, the spin length S, and a combined 

2 

quantity E at the quadrupole order. As definitions for 

2 2 

the quantities S and E, we use expressions which are 
completely analogous to the ones introduced at the pole- 
dipole order, see (jSTj) and (1531) . 



8 2 
— m 
ds 



6_ 2 

ds 



4'* 



S_ 2 

ds 



^ (u a S ab ^j + u b R ace b . d Q n dcae + 2 ?i aed u c + n ae u c u d ^J + R ace d 



I A ? 

d I z. abce ■ 

i - n UdUb 











ds ^ 


3 


I 






+ 2 








-( 


S 


ds \ 


3 






n l 


+ 2 






'i 


-( 


P 


ds \ 


3 








+ 2 





2 2 



s_ 

li 

8 2 S 2 



" ( cab £ 1 ^ £ ^j,2„ n n a f^bcde , o £ deftc 

ds V I ~ds 6 UaVbticde \n + 6n U 



3 



3^ 

2 2 , 16 2 / 
= 4 5a6P Q M 6 +y SabRcde a [i 



n bcde + n edc u b + 3 n bce u d + - n bc u d u e + - n de u b u c 



2 bcde + ledc u b + i n bce u d 



n 



— n de u b u c 
2 



1 2 
_ C 

2 



->u d u e + -n de u h u c \ 
2 2 / 



-| (^ a V b i? cde a + 2R cde a V b ip a ) (n bcde + n edc u b + 3 n bce u 



(112) 



(113) 



(114) 



(115) 



As becomes clear from (|112|) - (|115p these quantities are 
no longer conserved at the quadrupolar order. Note that 
this observation is independent of the choice of supple- 
mentary condition for the spin. A direct generalization 
of (|48|) - in terms of the quantities at the quadrupole 
order - only nullifies the first terms in (|112p - (|114l) . The 
conservation of the quantities in (|112[) - ([115[) depends on 
the details of the extended test body under considera- 
tion. Since we did not introduce any specific supplemen- 
tary condition for the quadrupole components, the lack 
of conserved quantities at the current order of approxi- 
mation is not unexpected. 



B. Comparison to Dixon's scheme 

Our starting point are the equations of motion given 
by Dixon in (13.7) and (13.8) of j27|, see also (168) and 
(169), as well as (171) and (172) in (28|. According to 



Dixon the spin and the center-of-mass equations of mo- 
tion - up to the quadrupole order - can be written in the 
form 

^-S ab = 2p[ Q u b l + -R cd J a I b ^ cde , (116) 
ds 3 

^rPa = ^RabcdU b S cd -]:S7 a R bcde I bedc : (117) 
ds 2 6 

where I abcd has the following symmetries: 

jabcd j(ab)(cd) jcdab (118) 

j{abc)d q jabcd _|_ jbcad jcabd q ^^9) 

Note that I abcd has no orthogonality properties with re- 
spect to u a , except the ones that can be deduced from 
its symmetries. In Dixon's work the equations of motion 
in (|116[) and (|117[) are also often written in terms of a 
different set of moments, termed J by him. To allow for 
a direct comparison, we provide the explicit transforma- 
tion rules between the I- and J-moments in appendix 

m 
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Orthogonal decomposition of (|116[) leads to 

n a n b — ^ cd - - n [a n b ^ R ,, f T 9cde 
ds ~~ ^PfPg n cde 1- , 



(120) 



p a bV h = S ab u b - -u b R c J a I»^ de . (121) 
With the definition m := u a p a equation (I12ip turns into 



p a = mu a + S ab u b - -u b R c J a I^ cde . 
Further, equation (|117|) can be written as 

A 1 1 

-rV — -7;ttcdb U b - - V b K c de * 
ds 2 3 



(122) 



(123) 



Insertion of (I122[) into f|l 23[) and orthogonal decomposi- 
tion of (|123[) leads to separate equations of motion for u a 
and m. This is completely analogous to the pole-dipole 
case, c.f. equation (l4"5l) and (j4T>)) . 

The form of Dixon's equations of motion in ([120j) and 
(|123[) is the most suitable one for a comparison with 
our results from the previous section. We introduce the 
orthogonal decomposition of I abcd : with its symmetries 
from (|118p and (|119[) already implemented, 

jdcab = Qdcab + 2 Q d < a u b ) + 2Q ab{d U C) + Q dc U a U b 

+Q ab u d u c - 2u {d Q c){a u b ^ , (124) 

where 4 



With the help of (|124[) and the momentum in terms of 
the Q-moments, i.e. 



p b = mu b + S ba u a + R a J ( -Q ab - • 



(130) 



the equations of motion (|120l) and (|123p turn into 



p a cP b d 



6S 



cd 



ds 
S 



P [ gP h }RdcJ 9 

+2Q ed u f u c ] 



-Q edcf + 4Q edf u c 



(131) 



—p" = --R ace »u*S ac + V d R ac 
as 2 



-Q dc 
3^ 



+-Q ,fce « + -Q ae£ V + i<3 dc u a u e 
3 3 2 



(132) 



Comparison of (|13ip to our equation of motion for the 

2 

spin (IllOp shows that the Q- and the n-moments only 
differ by a factor of 2, i.e. 

jrfcafc = 2 ^ l dcab + 2 ^dc(a u b) +2 lab{d u c) 

+ n dc u a u b + n ab u d u c - 2u ( - d n c ^ a u b ^ , (133) 



Qabcd 


— Q(ab)(cd 


Qcdab 


(125) 


)(abc)d 


= & 


Qabcd _|_ Qbcad _j_ Qcab( 


z = 0(126) 


Qabc 


= Q( Qb ) c , 




(127) 


Q(abc) 


= o 


Qabc + Qbca + gcab = 


0, (128) 


Q ab 


= Q (Qb) , 




(129) 



and all Q ab -- are orthogonal to the four velocity u a . No- 
tice that the symmetries of the Q-moments are identical 

2 

to the symmetries of the n- moments, which suggests that 
they are proportional to each other. 

4 Note that a similar decomposition has also been introduced in 
[2l . Therein quantities analogous to the Q abcd , Q abc , and Q ab 
are called stress-, flow-, and mass-quadrupolc. In contrast to 
our decomposition with respect to u a , Ehlers and Rudolph use 



as well as S ab —S ab '■ Furthermore, comparison of our 
equation of motion for the center-of-mass, i.e. equation 

__ 2 2 

(fTTTj) . with (fT3"2"]) yields m =m, and p a =P a . 



C. Energy- momentum tensor 

Our findings in the previous section allow for a com- 
pact representation of the energy-momentum tensor. In 
terms of Dixon's moments the decomposition in (|62l) - (|6~4")) 
takes the form 

the vector p a to perform the orthogonal decomposition. Further- 
more, in 12911 the decomposition was applied to the J-moments 
instead of the /-moments, see also (I A5E - 



±_dcab 



= -I 



dcab 



2 n ab{d u c) ~ n ab u d u c + 2u {d n°^ a u b) , 



^cau 



-S c{a u b) + — (2 n abc + n ab u c - 2 n c ( a u b) ) + A [ab) u c , 
ds V / 



2 2, 

3 



_ n b)cde + _ u b) n edc + 2 n b)cd u e + g ^c^d 



SAW 
ds 



(134) 
(135) 
(136) 
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Simplification by means of, e.g., 




-J V c (s c ^%)) + I J V d V c (I dcah 5 {i) ) . (138) 

Notice that this expression is not in canonical form, how- 
ever, it seems to be the simplest representation of the 
energy-momentum tensor density expressed in terms of 
the Dixon-moments. All constraint-type relations have 
been implemented into the energy-momentum tensor. 
Thus, all consequences of ([1]) can equivalently be written 
in a compact way as the equations of motion (|116|) and 
(|117p . the symmetries of the quadrupole (|118l) and (|119l) . 
and the energy- momentum tensor in the form (|138p . 



VII. CONCLUSIONS 

In this work we explicitly derived the equations of mo- 
tion for extended test bodies in General Relativity with 
the help of Tulczyjew's multipolar approximation method 
up to the quadrupolar order. To our knowledge this is 
the first time the method by Tulczyjew has been used 
beyond the pole-dipole order. In our derivation we put 
special emphasis on a transparent notation, which allows 
for a direct identification of the contributions of the mul- 
tipole moments at different orders, and explicitly carried 
out the canonicalization process. We would like to stress 
again, that we did not make any assumption for a sup- 
plementary condition in our derivation of the equations 
of motion. 

Our results are of direct relevance for other pertur- 
bation methods in the context of the general relativistic 
problem of motion. In particular the equations of motion 
in (IllOp and (|111[) , or - alternatively - the ones given in 
(|116| and (|117l) . as well as the energy- momentum tensor 
in (|138p , are needed in approximation schemes which aim 
for a description of self-gravitating compact objects and 
the gravitational radiation emitted by these systems. 



A. Structure of the equations of motion 



theorem B. We introduced the classification of two dif- 
ferent types of equations, termed "constraint" and "evo- 
lution" . As becomes clear from table U the same type 
of pattern of equations repeats at each multipolar order. 
In particular, one does not expect more than two equa- 
tions of the evolution-type in the context of Tulczyjew's 
approximation scheme. This can be viewed as support 
of Dixon's result [13, US HH , who provides a solution 
of what he calls the variational equations of mechanics^ 

- i.e. the combination of equation (fTJ and ([2]) - to any 
order. 

Our results in (|110p and (|111[) also allow for a direct 
comparison to the quadrupolar equations of motion de- 
rived in [lj| . The results therein were obtained via a dif- 
ferent multipolar approximation scheme which goes back 
to Papapetrou b 17|] . Although the same information is en- 
coded in the system of equations - Tulczyjew's as well as 
Papapetrou's method share the same starting point, i.e. 
the "conservation" of energy in the form of ([T]), further- 
more in both methods the full moments are taken into 
account - the final representation of the equations of mo- 
tion is different. This difference can mainly be ascribed to 
the use of the orthogonal decomposition of the moments, 
which is an integral part in Tuczyjew's procedure, and is 
introduced at a very early stage to support the deriva- 
tion of the canonical form. In particular the recursive 
transfer of higher order moments to lower differential or- 
ders in the canonical form, see the final result (|75[) at 
the quadrupolar order in this respect, yields structurally 
different equations. While it is mainly a question of prac- 
ticability which system of equations of motion should be 
given preference, the main benefit of the method by Tul- 
czyjew is its intrinsic covariance as well as its systematic 

- albeit somewhat laborious - way to generate a hierar- 
chical set of equations. 



B. Open problems 

Albeit the results obtained in this work are complete, 
in the sense that they cover the quadrupolar order in the 
context of Tulczyjew's approximation scheme in the most 
general way, there remain several interesting open ques- 
tions to be addressed. Some of these open problems are 
not specific to Tulczyjew's scheme, and are also inherent 
to other multipolar approximation schemes. 



C. Supplementary conditions 

As we have seen in sections [V] and IVI1 additional con- 
ditions are needed, starting at the pole-dipole order, to 



Table U gives a compact overview of the different equa- 
tions, which were obtained from the conservation law (fTJ) 
and the decomposition @ via repeated application of 



5 This notion goes back to Mathisson i, see 1 for more historical 
details. 
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TABLE I: Structure of the equations of motion. 



Order Quantity 



Type Equation 



Single-pole 



2(8-) 



c Gnj 

c US) 
e JT5) 



Dipole 



(ab)c 
(ab) 



o ab fo ab ,o a ,r 

o a (h ab ,o a ,t ab 
t (o ab ,o a ,t ab ,o a 



C 
C 
E+C 

c 

E 



El 



E2} 



Quadrupole 



(abc)d 
(abc) 



k { - ab)c (n ahc ,n 

l(a&) ^abc^ab^ 

n ab (n abcd ,n abc ,n ab ,n ab ,n a 
n a (n abcd ,n abc ,n ab ,n ab ,n a ^ 
n (n abcd ,ri abc ,n ab ,h abc ,h ab ,n a 



c 


m 


c 


m 


c 


(H2> 


c 


m 


E+C 


$E<3> 


c 


(87} 


1 E 





Type: "C"=constraint, "E" =evolution 



close the systems of the equations of motion. The ne- 
cessity for a spin supplementary condition can be seen 
as connected to the fixation of a specific representa- 
tive worldline inside the object. This is well-known in 
the special-relativistic context, where a spin supplemen- 
tary condition selects a representative worldline and vice- 
versa, see, e.g., (32j . In General Relativity, however, 
things are more subtle, c.f. [2!| [33l - l38j in this respect. 
In particular, it is difficult to prove that a spin supple- 
mentary condition fixes a representative worldline in a 
unique way. This is of course directly related to the fact, 

2 2 

that the quantities P a and S are combinations of geo- 
metrical quantities as well as of multipolar moments from 
different orders. Further studies are needed in the con- 
text of Tulczyjew's scheme when it comes to the choice 
of suitable supplementary conditions at higher orders. 

However, from a more pragmatic point of view every 
supplementary condition, which leads to a closed system 
of equations, can be used. The important question for 
applications is, whether the equations describe the mo- 
tion in an accurate way. 

The results derived in this paper could be interesting 



when studying whether two sets of moments contained in 
the energy-momentum tensor on the same, or infinitcsi- 
mally close worldlines, are equivalent. In the latter case 
one first has to shift the distributional energy-momentum 
tensor from one of the worldlines to the other. If two sets 
of moments are equivalent, then the coefficients in the 
canonical form of the corresponding energy-momentum 
tensors must be the same (by virtue of theorem B). We 
will not work this out in detail here. However, an imme- 
diate consequence is that an infinitesimal change of the 
worldline directly translates into an infinitesimal change 

of the moment n a . This has been used in [l], [H,[l6[ to ar- 
gue that, by a suitable choice of the representative word- 
line, one can restrict to n a (r) = for all r. In consid- 
eration of (EJTj) a change of h a corresponds to a change 

2 

of S ab Ub and will thus have an impact on the supple- 
mentary condition fulfilled by the spin. In this way, an 
infinitesimal change of the representative worldline can 
be related to an infinitesimal change of the spin supple- 
mentary condition. 



D. Combined and conserved quantities 

While it is possible to generate a hierarchical set of 
equations of motion, c.f. table HI with the help of Tul- 
czyjew's method, it is non-trivial to devise "combined" 
quantities at higher multipolar orders. A good example is 
the quantity I abcd as introduced in equation (|133[) . There 
is no straightforward algorithm in Tulczyjew's scheme 
which tells us how to construct it. While one may con- 
sider the introduction of quantities like I abcd as a mere 
question of taste - after all the main benefit is a more 
compact form of the equations of motion - it appears 
to be desirable to have such quantities at one's disposal. 
In particular when it comes to the search for conserved 
objects. 

As we have shown in section IVI A[ several of the com- 
bined quantities are no longer conserved at the quadrupo- 
lar order. This is of course directly linked to the fact that 
we did not assume any supplementary condition in the 
course of our derivation. Even the direct generalization 
and use of the supplementary conditions from the pole- 
dipole order does not yield a set of conserved quantities at 
the quadrupolar order. Our results in (11 12|) — (|115[) make 
clear, that additional conditions for the quadrupole mo- 
ments and/or symmetries of the underlying spacetime are 
needed to obtain a set of conserved quantities. 



E. Further approximations 

Due to the complexity of the expressions at the 
quadrupole order, further approximations are needed 
for applications. Most interesting is the restriction 
to some kind of mass-quadrupole, i.e., neglecting the 
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flow- and stress-quadrupole 6 . Corrections from the 
mass-quadrupole are needed, e.g, for the contributions 
quadratic in spin to the post-Newtonian dynamics, see 
[39l - |42l ]. The approach to higher order post-Newtonian 
spin dynamics in [43l - l46l | can incorporate quadrupole cor- 
rections in a straightforward way. Further, it is pos- 
sible to derive a canonical formalism from the energy- 
momentum tensor given in the present paper by the pro- 
cedure outlined in [13, HU in certain cases, c.f., |40| . 
Canonical methods also proved to be very useful in the 
post-Minkowskian approximation of single-pole objects, 
see, e.g., [49j]. Our results may also be used as input for 
perturbation methods [lOl - flij which aim for a description 
of systems with high mass ratios. Such methods have al- 
ready been applied to single-pole as well as to pole-dipolc 
objects. 



F. Regularization 

A straightforward application of the results in the 
present paper to self-gravitating compact objects is quite 
subtle from a mathematical point of view. For short, a 
strict mathematical definition of the product of distribu- 
tions does not exist, but would be required due to the 
non-linearity of Einstein's field equations. Therefore a 
distributional energy-momentum tensor as a source of 
the gravitational field makes mathematically no sense 
in General Relativity. However, this problem can be 
overcome, as in quantum field theory, by a regulariza- 
tion and renormalization program. In particular, dimen- 
sional regularization [50l l51j is most useful for theories 
involving gauge freedoms, like General Relativity. Di- 
mensional regularization has been employed successfully 
in post-Newtonian calculations [52h54| | to a high order of 
non-linearity. Notice that all results in the present paper 
hold for an arbitrary dimension of spacetime. 
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Appendix A: From /-moments to J-moments 

There is a one-to-one transformation from the I- 
moments to the J-moments, which are used interchange- 
ably in Dixon's work [2?], HH . At the quadrupole order 
it is explicitly given by 

jabcd j-[a[cb]c£] ^j-a[c|6|d] y-&[c|a|c£]^ (Al) 



jabed jd(ab)c ja(dc)b 

3 3 



(A2) 



where J has the following properties: 

jabcd ^j[ab][cd] jedab (A3) 

j[abc\d q jabcd _j_ jbcad _j_ jcabd q (A4) 

Thus, J abcd has the same (algebraic) symmetries as the 
Riemann tensor. Equation (|124p now becomes 

jabcd _ Q[a[cb]d] _ 2- lt [ a Q f, l [ cd l — 2u lC Q d ^ ah ^ 

-3u^ a Q b ^ c u d l (A5) 



In terms of J the equations of motion turn into 

-^-S ab = 2pK b ^ - ^R [a cde J b]cde , 
as 3 



(A6) 



^~Pa = \RabcdU b S Cd +\v a R bcde J hcde , (A7) 

as 2 6 

and the energy-momentum tensor becomes 



»6. 



(4) 



~ J V d V c (,J d ^ c 5 {i) 



(A8) 



Appendix B: Combined quantities in terms of o 
moments 

In this appendix we provide a summary of combined 

2 

quantities at the quadrupole order in terms of the o mo- 
ments, as used in the original decomposition of t abcd in 

©■ 



Notice that there can be slightly different ways to define such 
mass-, flow-, and stress-quadrupoles. They only have to ap- 
proach the correct Newtonian limit. 



jdcab _ 2 ^dcab + 2 2 dc(a u b) + 3 gah^c) 

2u {d o c ><°u 6 >) , (Bl) 



+ o dc u a u b + o ab u d u 
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S 



ab 



+u c t c[a u b] +2l 



g t [ab]c 

ds 



2u c u d 



S t cd[a 
ds 



, (B2) 



P b =m U b + U A + Race d (\ O abCe U d 

as \S 
+4 l bae u c u d + | l aec p b d + 2 l ae u c p b )j , (B3) 

mJt + S ab u a u b - ^R abc d o cba u d 



+U e Uf 



5 i t S 2 2 , (5 / j 2 , 
— t ef - — — t ef + 2— I u — t sef 
ds ds 2 ds V 9 ds 



+2uU cd{e R gcd f) 



(B4) 



Appendix C: Conventions &: Symbols 

In table|TI]we provide a list of symbols used throughout 
the text. Our convention for the signature of spacetime 
is —2. The curvature tensor is defined by (f2~Tj) . or equiv- 
alently by 

p b p b -p b _i_ p ep b p ep b /pi \ 

f^acd • — dc ,a 1 da ,c ~y L dc ea ^ da ^ ec • y^ 1 ) 

The delta function — 5^(x a — Y a (s)) is normalized 
as J J( 4 ) = 1. 



TABLE II: Directory of symbols. 



Symbol 



Explanation 



Geometrical quantities 



9ab 

9 

Ftij k 
yra 

u a 
s 

5(4) 



Metric 

Determinant of the metric 
Riemannian curvature 

Worldline within the worldtube of the body 
Velocity along the worldline Y a of the body 
Proper time along the worldline 
Four-dimensional delta function 



Matter quantities 





m 

i i 

P a , S c 
i i 
m, m 

l 

E 



2 2 

P a , s c 

2 2 

m, m 



2 

E 



jabcd jabc 

m,p a ,S ab 
Q ab ,Q abc , 



Energy-momentum tensor 

General multipole moments in the expansion 

of the energy momentum density 

Parts of the orthogonal decomposition of the 

general single-pole moment 

Parts of the orthogonal decomposition of the 

general dipole moment 

Parts of the orthogonal decomposition of the 

general quadrupole moment 

Canonical multipole moments 

Parts of the orthogonal decomposition of the 

canonical quadrupole moment 

Mass at the single-pole order 

Momentum and spin at the pole-dipole order 

Masses at the dipole order 

Conserved quantity at the dipole order (if 

spacetime has symmetries) 

Momentum and spin at the quadrupole order 

Masses at the quadrupole order (not conserved 

in general) 

"Combined" quantity at the quadrupole order 

(not conserved in general) 

Dixon's quadrupolar moments 

Dixon's "combined" quantities 

Parts of the orthogonal decomposition of 

Dixon's quadrupole moment I" ° 



Operators 



Vi, 
5/ds 



Covariant derivative 
Total covariant derivative 
Partial derivative 
Spatial projector 



Accents 



Denotes the density of an object 

Denotes the orthogonal projection of an index 

Denotes the derivative 5/ds 
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